In [1]:

    

library("NADA")


    

    




Loading required package: survival

Attaching package: 'NADA'

The following object is masked from 'package:stats':

    cor

In [2]:

    

main_columns = c(
    'epa_rain_zone', 'season', 'primary_landuse',
    'parameter', 'fraction', 'units', 'res', 'qual'
)
df = read.csv("nsqdata_example_subset.csv")[main_columns]
head(df)


    

    
Out[2]:





epa_rain_zone
season
primary_landuse
parameter
fraction
units
res
qual

	
1
1
spring
Residential
Cadmium
Total
ug/L
10
<
	
2
1
spring
Residential
Cadmium
Total
ug/L
10
<
	
3
1
spring
Residential
Cadmium
Total
ug/L
10
<
	
4
1
spring
Residential
Cadmium
Total
ug/L
10
<
	
5
1
spring
Residential
Cadmium
Total
ug/L
1
<
	
6
1
spring
Residential
Cadmium
Total
ug/L
1
<

In [3]:

    

dfros = NADA::ros(df$res, df$qual=='<')
median(dfros)


    

    




Warning message:
In NADA::ros(df$res, df$qual == "<"): Dropped censored values that exceed max of uncensored values.





    
Out[3]:




0.351893483458947